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Systems  of  Evolution  Equations  in  Thermomechanical  Processes 
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Research  Description:  In  the  framework  of  continuum  mechanics  the  motion  of  de¬ 
forming  bodies  is  described  by  systems  of  nonlinear  evolution  equations,  that  arise  by 
combining  balance  laws  with  constitutive  relations  characterizing  the  material  response. 
The  nonlinear  character  of  the  material  response  often  induces  a  destabilizing  mechanism 
that  competes  with  dissipative  mechanisms,  such  as  viscosity  or  thermal  diffusion.  As  a 
result  of  the  competition  coherent  structures  may  appear,  which  at  some  level  of  modeling 
manifest  themselves  as  singularities  in  the  solutions  of  the  corresponding  model.  These 
structures  are  diverse  in  nature,  ranging  from  shock  waves  to  shear  bands  to  propagating 
phase  boundaries,  and  no  common  theory  can  encompass  all  of  them  at  present.  Various 
instances  of  such  phenomena  were  studied  a^  part  of  this  project: 

Shear  Bands.  One  of  the  most  striking  manifestations  of  instability  in  solid  mechanics  is 
the  localization  of  shear  strain  into  narrow  bands  during  high  speed,  plastic  deformations 
of  metals.  Shear  bands  are  structures  with  a  sharp  transition  layer  in  the  displacement 
field  and  a  corresponding  concentration  in  strain.  As  they  often  precede  rupture  their 
understanding  is  critical  for  the  improvement  of  materials,  and  their  study  has  attracted 
a  lot  of  attention  from  diverse  points  of  view.  The  objective  of  this  program  is  to  under¬ 
stand  a  thermoplastic  instability  mechanism  proposed  for  the  explanation  of  shear  band 
formation  at  high  strain  rates.  According  to  this  scenario,  two  main  factors  are  competing 
during  shear  band  initiation:  First,  effective  strain-softening  response,  resulting  eis  the 
net  outcome  of  the  influence  of  thermal  softening  on  the  (under  isothermal  conditions) 
strain-hardening  response  of  metals.  This  is  a  destabilizing  factor  having  the  tendency 
to  amplify  nonhomogenuities.  And  second,  strain-rate  hardening  which  heis  the  opposite 
effect  of  diffusing  nonhomogenuities.  To  understand  this  competition  various  test  defor¬ 
mations  were  considered.  They  lead  to  systems  of  nonlinear  partial  differential  equations. 
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which  axe  studied  by  means  of  nonlinear  analysis  techniques.  The  techniques  include 
the  use  of  comparison  estimates,  the  existence  of  invariant  regions  and  exploiting  scale- 
invariance  properties  that  specific  models  are  equipped  with.  Some  of  the  main  findings 
are  summarized  below  (for  details  see  the  technical  description  part  of  the  report). 

To  assess  the  contributions  of  thermcd  softening,  strain  hardening  and  strain-rate  sensitiv¬ 
ity,  a  shearing  deformation  caused  by  prescribed  tractions  has  been  analyzed  for  a  model 
based  on  a  power  law  and  for  which  explicit  thermal  effects  are  taken  into  account.  It  turns 
out  that  the  parameter  space  formed  by  the  powers  can  be  decomposed  into  three  distinct 
regions  across  which  the  behavior  of  the  corresponding  solutions  changes  drastically;  from 
blow  up,  to  weak  instability,  to  asymptotic  stability.  The  first  two  regions  are  associated 
with  the  development  and  evolution  of  spatial  nonhomogenuities  in  the  strain,  in  the  form 
of  shear  bands,  and  are  characterized  by  a  collapse  of  the  ability  of  the  material  to  diffuse 
the  applied  stresses  [1,  3]. 

To  understand  the  mechanics  of  development  and  propagation  of  shear  bands,  a  shearing 
deformation,  caused  by  steady  boundary  shearing,  of  a  strain-softening,  rate-dependent 
material  has  been  analyzed.  A  linearized  stability  analysis  of  the  time-dependent  lin¬ 
earized  equations  characterizes  the  necessary  degree  of  strain  softening  to  produce  growing 
perturbations  of  the  uniform  shearing  solutions.  In  a  result  for  the  nonlinear  problem,  the 
evolution  of  an  initial  large  concentration  of  strain  is  monitored  and  shown  to  produce  a 
fully  developed  shear  band  at  large  times;  outside  the  band  the  material  unloads  [2,  4]. 


Self-Similar  Limits.  The  invariance  of  hyperbolic  systems  of  conservation  laws  under 
dilations  of  coordinates  is  a  key  property  underlying  much  of  the  current  theory.  Viscous 
perturbations  introduce  an  additional  parabolic  scale,  and  this  invariance  is  lost.  At  present 
it  is  not  known  how  these  two  scales  interact  in  the  context  of  viscous  limits  for  the  Riemann 
problem.  One  illuminating  step  in  this  direction  is  to  study  artificial  regularizations,.' 
rigged  so  as  to  preserve  the  invariance  under  dilatations  of  coordinates.  It  leads  to  the 
study  of  singular  perturbations  for  non- autonomous  boundary-value  problems,  and  the 
limiting  process  involves  variation  estimates.  This  approach  was  tested  for  the  equations 
of  isentropic  gas  dynamics  in  Eulerian  coordinates  [7],  and  is  currently  under  study  for 
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larger  systems.  Recently,  a  similar  idea  was  employed  for  the  case  of  the  Broadwell  system  Codes 
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with  Riemann  data.  The  system  is  modified  so  as  to  incorporate  the  invariance  under 
dilations  of  the  limiting  conservation  laws,  the  resulting  self-similar  solutions  are  studied, 
and  the  fiuid-dynamic  limit  is  carried  out  [8]. 

Viscoelasticity.  A  simple  model  problem  was  studied  intended  to  test  the  following 
conjecture:  That  lack  of  monotonicity  in  the  steady  shear  stress  versus  steady  shear  strain- 
rate  leads  to  non-Newtonian  response,  manifested  by  the  occurrence  of  spurt  in  viscometric 
fiows  of  certain  pclymer  solutions.  It  was  shown  that  a  class  of  discontinuous  steady  states 
(associated  with  spurt)  are  asymptotically  stable  [5]. 

Technical  Description  A  technical  description  of  the  projects  supported  by  the  U.S. 
Army  Research  Office  under  Grant  DAAL03-88-K-01S5  follows.  We  include  studies  that 
are  in  the  preparation  process  and  which  (when  published)  will  acknowledge  support  from 
this  Grant. 

Studies  of  the  phenomenon  of  shear  band  formation  at  high  strain  rates 

[1]  Strain  softening  in  viscoelasticity  of  the  rate  type  J.  Integral  Equations  Appl.  3  (1991), 
195-238. 

In  this  article  we  study  the  behavior  of  solutions  of  the  system  of  partial  differential 
equations 

Vt=<Tx 

(1) 

Ut  =  Vx 

(2)  <7  =  t(u)v^ 

where  t(u)  satisfies  r(u)  >  0  and  t'(u)  <  0.  The  equations  describe  the  plastic  shearing  of 
a  plate  of  unit  thickness.  In  this  context,  v  stands  for  the  shearing  velocity,  u  for  the  plastic 
strain  and  a  for  the  shear  stress.  Our  objective  was  to  use  this  model  as  a  paradigm,  in 
order  to  analyze  the  competition  between  the  destabilizing  effect  of  strain  softening  versus 
the  stabilizing  influence  of  strain-rate  sensitivity,  the  strength  of  which  is  measured  by  the 
(typically  small)  parameter  n. 
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A  general  existence  and  continuation  theory  of  classical  solutions,  for  problems  fitting 
under  the  structure  (1—2),  is  carried  out.  Subsequently,  the  time-evolution  of  the  solutions 
corresponding  to  two  types  of  loading  is  studied: 

Prescribed  tractions  ^(O,^)  =  <t(1,<)  =  1:  A  characterization  of  global  existence  or  blow 
up  of  solutions  in  terms  of  convergence  properties  of  the  integral  is  first  es¬ 

tablished.  Then  the  response  of  solutions  is  discussed  under  various  hypotheses  for  the 
function  r(u).  We  present  the  result  for  a  power  law  r(u)  =  l/u"*  with  m  >  0.  The 
parameter  region  can  be  decomposed  into  three  distinct  subregions  accross  which  the  re¬ 
sponse  changes  drastically.  For  0  <  m/n  <  1/2,  solutions  are  asymptotically  attracted  to 
a  stable  time-dependent  state;  forl/2<m/n<l,  this  state  loses  its  attracting  property; 
and  for  m/n  >  1,  solutions  blow  up  in  finite  time.  This  change  of  behavior  is  associated 
with  inability  of  the  material  to  diffuse  the  applied  stresses,  and  developement  of  spatial 
nonuniformities  in  the  strain  (associated  with  shear  band  formation)  at  the  boundaries. 
Prescribed  velocities  u(0,t)  =  0,  v(l,t)  —  1:  We  consider  a  power  law,  r(u)  =  l/u'",  and 
use  scaling  invariance  properties  together  with  the  existence  of  invariant  regions  for  an 
associated  parabolic  system  to  discuss  the  stability  of  a  special  class  of  solutions,  describing 
shearing.  It  turns  out  that,  for  0  <  m/n  <  1,  the  uniform  shearing  solutions  are  stable 
and  every  solution  is  eisymptotically  attracted  as  t  >  oo  to  one  of  them. 


[2]  Shear  strain  localization  in  plastic  deformations  in  "Shock  Induced  Transitions  and  Phase 
Structures  in  General  Media”,  (R.  Fosdick,  E.  Dunn  and  M.  Slemrod,  eds),  IMA  Volumes 
in  Mathematics  and  its  Applications,  Springer  Verlag,  1991  (to  appear). 


In  this  article  the  system  (1  —  2)  above  is  discussed.  The  first  part  is  a  summary  of  the 
results  presented  in  [1]. 

In  the  second  part,  we  take  up  a  power  law  a  =  (l/u’")t;”  for  shearing  caused  by  prescribed 
velocities  and  consider  the  parameter  region  m/n  >  I  (complementary  to  the  Ccise  discussed 
in  [1]).  It  is  shown  that,  for  initial  data  corresponding  to  a  linear  velocity  profile  and  a 
concentration  of  strain,  either  the  corresponding  strain  blows  up,  or  as  t  oo 


0{t-^) 
l  +  0(t-*) 


X  €  [0, 1/  -  <5] 

a:  €  [y  +  ^,  1] 
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and  u(x,t)  converges  outside  the  region  [t/  —  S,  t/  +  SJ  to  a  limiting  configuration.  Thus  the 
material  unloads  outside  the  band. 

[3]  On  adiabatic  shear  bands  (in  preparation). 

This  manuscript  is  still  in  the  preparation  process.  Here  we  present  some  results  available 
at  the  time  of  the  writing. 

Our  goal  in  this  article  is  to  extend  the  understanding  of  the  shear  band  formation  problem 
obtained  in  [1,2]  for  cases  that  explicit  thermal  effects  are  taken  into  account.  We  consider 


the  model 

ci 

II 

Ci 

(3) 

Ut  =  Vx 

9t  =  avx 

(4) 

a  =  fi{9,u) 

with  boundary  conditions:  either  prescribed  tractions  or  prescribed  velocities.  The  func¬ 
tion  jj.  is  assumed  decreasing  in  9^  and  n  is  positive  measuring  the  degree  of  strain-rate 
sensitivity.  Regarding  the  dependence  of  //  in  u  both  strain  hardening  (^u  >  0)  and  strain 
softening  <  0)  are  admitted  as  possibilities,  but  an  assumption  that  guarantees  that 
thermal  softening  dominates  over  strain  hardening  is  instead  imposed.  This  resulting  sys¬ 
tem  is  a  model  to  assess  the  interplay  of  thermal  softening,  strain  hardening  and  strain-rate 
sensitivity  on  the  response  of  pleistic  flows. 

For  loading  effected  by  prescribed  tractions,  a  characterization  of  global  existence  or  blow¬ 
up  in  terms  of  the  behavior  of  lJ-{9,u)  is  presented.  Then  the  behavior  of  solutions  is 
analyzed.  To  illustrate,  consider  the  case  of  a  power  law  a  =  9''u'^u^,  with  parameters 
1/  <  0,  m  and  n  >  0.  It  turns  out  that  the  parameter  region  can  be  decomposed  into 
three  subregions:  (a)  i/  -f  m  <  — n,  (b)  — n  <  u  +  m  <  — n/2  and  (c)  — n/2  <  u  +  m  <  0. 
In  region  (c)  the  solutions  exhibit  a  stable  response,  in  region  (b)  they  develop  spatial 
nonuniformities  and  in  region  (c)  they  blow  up.  This  generalizes  the  analysis  in  [Ij. 

In  the  case  of  prescribed  velocities,  the  case  when  strain  dependence  is  absent  is  considered 
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e,  =  (-(«)<+■ 

where  fi{d)  is  decreasing  to  0.  It  is  shown  that,  for  certain  classes  of  initial  data,  either 
the  solution  blows  up  or  the  uniform  shearing  solution  is  unstable. 

i 

[4]  Nonlinear  analysis  techniques  for  shear  band  formation  at  high  strain  rates,  Appl.  Mech. 
Reviews  (3)  45  (1992),  82-94. 

This  article  presents  a  review  of  the  results  of  [1-3]  emphasizing  the  ramifications  on  the 
mechanics  aspects  of  the  shear  band  formation  problem.  In  addition,  a  linearized  stability 
analysis  for  the  uniform  shearing  solutions  in  the  context  of  (1  —  2),  taking  account  of  the 
time  dependence  of  the  uniform  shearing  solutions,  is  presented  here. 

Studies  related  to  viscoelasticity 

[5]  Stability  of  discontinuous  steady  states  in  shearing  motions  of  a  non- Newtonian  fluid  (with 
J.A.  Nohel  and  R.L.  Pego),  Proc.  Roy.  Soc.  Edinburgh  115A  (1990),  39-59. 

We  study  the  nonlinear  stability  of  discontinuous  steady  states  of  the  following  model 
initial-boundary  value  problem  in  one-space  dimension,  describing  incompressible,  isother¬ 
mal  shear  flow  driven  by  a  constant  pressure  gradient 

=  5x  =  (O-  -f  Uj  -I-  fx)x 
<rt  +  (T  =  g(vj;) 

with  boundary  conditions  5'(0,  t)  =  0 ,  u(l,<)  =  1.  The  non-Newtonian  contribution  to  the 
stress  satisfies  a  differential  constitutive  law.  The  problem  shares  certain  key  features  with 
models  proposed  to  explain  the  occurrence  of  ’’spurt”  phenomena  in  non- Newtonian  flows. 
This  system  admits  steady  state  solutions  (u(x),a(z))  satisfying 

ff(vx)  +  Vx  +  f^  =  0  0- =  g(vx) 

In  the  ca.se  of  interest  the  function  w{^);  =  g{^)  -f-  Os  non-monotone,  and  as  a  result  the 
problem  admits  multiple  steady  states,  leading  to  velocity  profiles  with  kinks.  We  show 
that  every  solution  tends  to  a  steady  state  as  t  —*  oo.  Moreover,  we  show  that  the  steady 


6 


states  for  which  Vx(x)  takes  values  in  the  monotone  increasing  parts  of  the  graph  of 
aie  stable. 


[6]  Nonlinear  stability  in  non-Newtonian  flows,  (with  J.A.  Nohel  and  R.L.  Pego),  in  ’’Mul¬ 
tidimensional  Hyperbolic  Problems  and  Computations”  (J.  Glimm  and  A.  Majda  eds.), 
pp.  251-260,  IMA  Volumes  in  Mathematics  and  its  Applications,  Vol.  29,  Springer  Verlag, 
New  York,  1991. 

This  is  a  conference  proceedings  report  outlining  the  results  in  [6]. 


Studies  related  to  self-similar  limits 

[7]  A  limiting  viscosity  approach  for  the  Riemann  problem  in  isentropic  gas  dynamics,  (with 
M.  Slemrod),  Indiana  Univ.  Math.  J.  38  (1989),  1047-1074. 

We  consider  the  equations  of  isentropic  gas  dynamics  in  Eulerian  coordinates 

Pt  +  {pu)x  =  0 

(6) 

(p«)«  +  +p(p))*  =  0 

with  p'{p)  >  0  and  Riemann  data.  We  construct  solutions  of  this  problem  as  limits  of 
solutions  of  a  ’’viscosity”  regularized  problem  that  is  rigged  so  as  to  preserve  the  invari¬ 
ance  of  (6)  under  dilations  of  the  independent  variables.  The  limiting  procedure  is  based 
on  variation  estimates  for  the  solutions  of  an  associated  system  of  ordinary  differential 
equations.  The  solutions  thus  constructed  may  contain  vacuum  regions.- 

[S]  Self-similar  fluid  dynamic  limits  for  a  modified  Broadwell  system  (with  M.  Slemrod), 
(preprint). 


We  study  self-similar  hydrodynamic  limits  for  a  modified  Broadwell  model,  which  once 
transformed  to  self-similar  variables  reads: 

fi{±oo)  =  fi,±,i  =  1,2,3 
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with  data  corresponding  to  ’’Maxwellians”,  that  is  -  fi,±f2,±  =  0.  The  variation 
estimates  as  £  \  0  are  beised  on  the  observation  that  Q  =  —  fif2  can  only  vanish  at  the 

singular  points  ^  =  — 1,0, +1  (corresponding  to  characteristic  directions  for  the  Broadwell 
model).  Thus  the  /i’s  can  change  monotonicity  only  at  these  points.  This  is  the  main 
observation  that  allows  the  hydrodynamic  limit  £  \  0  to  be  carried  out  for  any  data. 
Because  of  the  singulcirities  at  — 1,+1  and  0,  the  existence  of  solutions  for  (B)  requires 
some  work.  The  approach  used  is  based  on  linearization,  uses  the  Frobenius  theory  for 
systems,  and  is  restricted  for  data  that  are  close  to  a  Maxwellian. 
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